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A discussion of results for short and long-range multiplicity correlations (forward-backward) are 
presented for several systems (AuH-Au, CuH-Cu, and pp) and energies (e.g. ^snn = 200, 62.4, and 
« 20 GeV). These correlations are measured with increasing values of a gap in pseudorapidity, 
from no gap at midrapidity to a separation of 1.6 units (JT] | = 0.8). For the highest energy, central 
Ah-A collisions, the forward-backward correlation strength maintains a constant value across the 
measurement region. In peripheral collisions, at lower energies, and in pp data, the maximum 
appears at midrapidity. This result may indicate the possible formation of high density matter for 
central Ah-A collisions at ^/snn = 200 GeV. 
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1. Introduction 

The study of long-range multiplicity correlations has a long history in e+ — e^ and p — p{p) 
experiments [|TJ, ^. More recently, these studies have been expanded to the realm of ultra-relativistic 
heavy ion collisions at RHIC [^. Correlations that extend over a large range in pseudorapidity (t]) 
have the potential to probe the early stages of heavy ion collisions. In hadron-hadron experiments, a 
Unear relationship relating the multiplicity in the forward direction (Nf) to the average multiplicity 
in the backward direction (Nb) was found [Q]: 

<Nt{Nf)>=a + bNf (1.1) 

The linear coefficient, b, is referred to as the correlation coefficient (or forward-backward cor- 
relation strength). The value, a, is a measure of uncorrected particles. For the definition provided 



in Eq. ( |1.1| ) the maximum value of the forward-backward correlation strength, b, is one. This 
corresponds to complete correlation between the produced particles. For b = 0, all particles are un- 
corrected. The relationship from Eq. ( |1.1[ ) can be expressed in terms of the following expectation 
values [0]: 

^_ <NfNh>-<Nf><Nb> _Dlf 
<NJ>-<Nf >2 Dj^ 

where D| r and D^^ are the backward-forward and forward-forward dispersions, respectively. This 
result is exact and model independent [^. All quantities in Eq. ( [T!^ ) are directly measurable. The 
large pseudorapidity and angular coverage of the STAR experiment at RHIC provides an ideal 
means to measure particle correlations with precision [^. 

In phenomenological particle production models, the hadronization of color strings formed by 
partons in the collision account for the produced soft particles. In the Dual Parton Model (DPM), 
at low energies, these strings are formed between valence quarks ^. At higher energies, there are 
additional strings formed due to the increasing contribution from sea quarks and gluons. In the 
DPM, this production mechanism is predicted to produce long-range multiplicity correlations in 
pseudorapidity. A long-range correlation could be an indication of multiple elementary inelastic 
collisions in the initial stages of an A-i-A collision [^. Additionally, the Color Glass Condensate 
(CGC) model predicts long-range correlations in association with the formation of an intermediate 
state, the Glasma [^. The CGC model is a description of high density gluonic matter at small-x. In 
the CGC model, the incoming Lorentz contracted nuclei are envisioned as two sheets of colliding 
colored glass |^. In the initial stage of the collision, the color electric and magnetic fields are 
oriented transverse to the beam direction (and each other). After the collision there are additional 
sources formed that produce longitudinal color electric and magnetic fields. This state is referred 
to as the Glasma, an intermediate phase between the CGC and Quark Gluon Plasma (QGP) [p^]. 
The longitudinal fields provide the origin of the long-range correlation and are similar to the strings 
as formulated in the DPM. 

2. Data Analysis 

The data utilized for this analysis was acquired by the STAR (Solenoidal Tracker at RHIC) 
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Figure 1: Schematic diagram of the measurement of a forward-backward correlation. 



experiment at the Relativistic Heavy Ion Collider (RHIC) [Q]. RHIC has delivered collisions of the 
heavy nuclei gold (Au) and copper (Cu), as well as proton-proton collisions. The data discussed 
here is from Au-i-Au, Cu-i-Cu, and pp collisions. For Au-i-Au collisions, the center-of-mass colli- 
sion energies i^/sNN) used are 200 and 62.4 GeV; Cu-i-Cu, ^/s^n = 200, 62.4, and 22.4 GeV; and 
pp, ^ysNN = 400, 200, and 62.4 GeV. The main tracking detector at STAR is the Time Projection 



Chamber (TPC) [11 ]. The TPC is located inside a solenoidal magnet generating a constant, longi- 
tudinal magnetic field. For this analysis, data was acquired at the maximum field strength of 0.5 T. 
All charged particles in the TPC pseudorapidity range 0.0 < |t] | < 1.0 and with pr > 0.15 GeV/c 
were considered. This rj range was subdivided into forward and backward measurement intervals 
of total width O.lrj. The forward-backward correlations were measured symmetrically about T] = 
with varying rj gaps (measured from the center of each bin), At] = 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 
1.6, and 1.8 units in pseudorapidity. This is different than the standard two-particle correlation 
function, which utihzes a relative Tj difference for particle pairs. Figure [l] depicts a schematic rep- 
resentation of the forward-backward correlation and its measurement. The collision events were 
part of the minimum bias dataset. The minimum bias collision centrality was determined by an of- 
fline cut on the TPC charged particle multiplicity within the range JT] | < 0.5. The centralities used 
in this analysis account for 0-10, 10-20, 20-30, 30-40, 40-50, 50-60, 60-70, and 70-80% of the total 
hadronic cross section. An additional offline cut on the longitudinal position of the collision vertex 
(Vj) restricted it to within ±30 cm from z = (center of the TPC). Approximately two million 
minimum bias events at ^/snn - 200 GeV and 0.7 million events at ^/s^n = 62.4 GeV satisfied 
these requirements and were used for this analysis. Corrections for detector geometric acceptance 
and tracking efficiency were carried out using a Monte Carlo event generator and propagating the 
simulated particles through a GEANT representation of the STAR detector geometry. 

In order to eliminate the effect of statistical impact parameter (centrality) fluctuations on this 
measurement, each relevant quantity i{Nf), {Nh), {Nf) , and (^NfNb)) was obtained on an event- 
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Figure 2: Forward-backward correlation strength, b, as a function of the pseudorapidity gap (Arj) in 0-10% (black 
squares) and 40-50% (blue squares) most central ^jjvw = 200 GeV Au+Au data. If only short-range correlations were 
present, the expectation is a rapidly falling b value with increasing At] , as seen in the 40-50% most central data. The 
approximately flat behavior of b in 0-10% data indicates the presence of a long-range correlation. 

by-event basis as a function of STAR reference multiplicity, Ndi- A linear fit to {Nf) and {Nb), 
or a second order polynomial fit to {Nf) and (^NfNt) was used to extract these quantities as 
functions of Nch- Tracking efficiency and acceptance corrections were applied to each event. Due 
to statistical limitations, it is not possible to apply corrections for every value of Nch- One value for 
the correction, calculated for each centrality, is applied to the resultant values (A^/) , {Nb) , {N/) , 
and (^NfNb) for every event in that centrality. Therefore, all events falling within a particular 
centrality have the same correction. These were then used to calculate the backward-forward and 
forward-forward dispersions, Dj,- andD^ ,, binned according to the STAR centrality definitions and 
normalized by the total number of events in each bin. A strong bias exists if the centrality definition 
overlaps with the measurement interval in T] space. Therefore, there are three T] regions used to 
determine reference multiplicity: |t]| < 0.5 (for measurements of At] > 1.0), 0.5 < |t]| < 1.0 (for 
measurements of At] < 1.0), or |t]| < 0.3 + 0.6 < |t]| < 0.8 (for measurement of At] = 1.0). This is 
possible due to the flat T] yield of the STAR TPC within the measurement ranges being considered. 



3. Results 

Figure § shows the results for the correlation strength b for most central (0-10%) and mid- 
peripheral (40-50%) Au-i-Au collisions at ^/sm = 200 GeV. The short-range correlation (SRC) 
is defined as the correlation for particle separations less than At] = 1 .0. Long-range correlations 
(LRC) are considered beyond Arj > 1 .0. Figure ^ shows that central ^/s]^ = 200 GeV Au-i-Au data 
exhibits a strong LRC. The correlation is approximately flat as a function of the pseudorapidity 
gap. More peripheral (40-50%) Au-i-Au data shows only a SRC that decreases rapidly as a function 
of the pseudorapidity separation. It is expected that if only short-range correlations are present 
(mostly from sources such as cluster formation, jets, or resonance decay), then b would decrease 
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Figure 3: The centrality dependence of the FB correlation strength for Au+Au at ^/Inn = 200 GeV. Most central (0- 
10%) is black; 10-20%, red; 20-30%, green; 30-40%, magenta; and 40-50%, blue. The long-range correlations persist 
for the most central Au-l-Au data, but is similar to the expectation of only short-range correlations at 30-40% most central 
data. 



rapidly as a function of Arj . It has been argued that the presence of these long-range correlations at 
larger than Arj = 1.0 cannot be due to final state effects, and are therefore due to the conditions in 
the initial state ^. This may be phenomenologically similar to fluctuations larger than the event 
horizon being created during an inflationary phase in the early universe that later become smaller 



than the event horizon 1 12]. 



The evolution of the forward-backward (FB) correlation strength as a function of centrality 
in ^/snn = 200 GeV Au-i-Au data is presented in Figure ^. There is a smooth evolution in the 
trend from most central (0-10%) to mid-peripheral (40-50%) that shows a decrease of the FB 
correlation strength as a function of collision centrality. Were only short-range correlations present, 
the behavior for all centralities should resemble the 40-50% result. There appear to be long-range 
correlations for the top 30% central Au-i-Au data at ^/snn = 200 GeV. 

The system size dependence of the FB correlation strength has been studied in Au-i-Au, CuH-Cu, 
and pp collisions. Figure ^ shows the comparison of the FB correlation strength for 0-10% most 
central ^/s^n = 200 GeV Au-i-Au and Cu-i-Cu. Central ^smn = 200 GeV Cu-i-Cu exhibits littie 
difference compared to Au-i-Au. Central Cu-i-Cu still exhibits a long-range correlation. 

Figure ^ shows the system size dependence of the FB correlation strength for mid-peripheral 
(40-50%) Au-i-Au, Cu-i-Cu, and pp data at ^/sj^ = 200 GeV As a function of the T] gap both 
the heavy ion and pp result exhibit a similar evolution. There is good agreement between all three 
systems within the calculated errors. The data point closest to mid-rapidity is predominantly driven 
by short-range correlations. This is exhibited in the ordering of the data near mid-rapidity. The FB 
correlation strength for AuH-Au is largest, followed by CuH-Cu, and pp, following a similar trend in 
multiplicity. 

In addition to the system size dependence of the FB correlation strength, the energy depen- 
dence has also been studied. Data is available at several different energies for both heavy ions 
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Figure 4: A comparison of the FB correlation strength in 0-10% most central Au+Au (squares) and Cu+Cu (triangles) 
data at ^/s/^ = 200 GeV. The long-range correlation seen in Figure g persists at a slightly lower value in Cu+Cu at the 
same energy. 
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Figure 5: The system size dependence of the FB correlation strength for 40-50% most central Au+Au (squares) and 
Cu+Cu (triangles) at y'SNN = 200 GeV. Also shown are the results for minimum bias pp (red circles) at the same energy. 
The mid-peripheral results for all three systems show a similar evolution as a function of the T] gap. 



(Au+Au and Cu-i-Cu) and pp. The energies available for the analysis of the FB correlation strength 
in pp collisions are ^/snn = 400, 200, and 62.4 GeV. For Au+Au and Cu-i-Cu the available energies 
are ^/s^n = 200, 62.4, and « 20 GeV. Figure ^ is the energy dependence of the FB correlation 
strength for all three energies of pp. There is good agreement between the FB correlation strength 
at y/sf^ = 400 and 200 GeV. The ^/s]^ = 400 GeV data shows a larger short-range correlation, but 
plateaus at the same value (b = 0.05) as the ^/sNN = 200 GeV result for Arj > 1.0. In contrast, it 
is seen that at y'siviv = 62.4 GeV the FB correlation strength is smaller across the entire T] region. 
Unlike the other two energies, the FB correlation strength at this lower energy approaches b = at 
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Figure 6: Energy dependence of the FB correlation strength for pp data, as a function of Arj. Three energies are 
shown: y'Siv/v = 400 (green circles), 200 (black circles), and 62.4 GeV (red circles). The y'S^w = 400 and 200 GeV 
results exhibit close agreement for large values of At]. The ^/snn = 62.4 GeV FB correlation strength is lower than the 
other two and goes smoothly lo h = 0. 



large values of Arj . 

The energy dependence of Au+Au and Cu+Cu will be discussed briefly. In the most central 
data (0-10%) for both systems at ,,/snn = 200 and 62.4 GeV, there is a plateau of the FB correlation 
strength for Arj > 1.0. However, the plateau at ^/snn = 62.4 GeV is w 30% lower in Au+Au and w 
40% lower in Cu+Cu. The FB correlation strength for Cu+Cu at ^/sNN = 22.4 GeV is quantitative 
agreement w/ the result at ^/s^n = 62.4 GeV. The FB correlation for 40-50% most central data at 
all energies for both systems is in qualitative agreement. They all exhibit behavior consistent with 
only short-range correlations. 

4. Discussion 



The Dual Parton Model (DPM) ascribes the long-range correlation to fluctuations in the num- 
ber of elementary, inelastic collisions, dictated by unitarity [^. It was shown that the FB correlation 
strength for pp and mid-peripheral (40-50%) heavy ion data is dominated by short-range correla- 
tions. This contrasts with more central heavy ion data that exhibits a large value of the FB correla- 
tion strength for Arj > 1.0. The experimental measurement of a long-range correlation is indicative 
of multiple partonic interaction in the case of central and semi-central heavy ion interactions. 

One possible explanation of the centrality dependence of the FB correlation strength as shown 
in Figure || relates to the idea of the formation of a central, partonic core surrounded by a hadronic 
corona. This has been discussed in terms of hydrodynamic properties and finite formation time 



of the strongly interacting quark gluon plasma (sQGP) [13, 14], as well as attempts to quantify 



individual contributions from the core and corona [15]. If the long-range correlation manifests 
due to partonic effects (as described in the DPM), the partonic core formed in heavy ion collisions 
would be the source of this correlation. It is expected from geometrical arguments that this partonic 
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core would have a larger volume in central heavy ion collisions, compared to more peripheral 
overlap. Therefore, the increasing influence of a partonic corona would lead to a larger long- 
range correlation. The evolution with centrality also agrees with predictions from the DPM and 
CGC/Glasma phenomenologies. ^ 

5. Summary 

The energy and system size dependence of the FB correlation strength, b, has been discussed. 
It has been studied for three colliding systems: Au+Au, Cu+Cu, and pp. The pp energies range 
from y/SNN = 400 to 62.4 GeV, while those for heavy ions range from ^/snn = 200 to ss 20 GeV. For 
the most central heavy ion collisions at ^/siviv = 200 GeV, the value of the FB correlation strength 
is approximately flat across a wide range in Arj . This behavior slowly disappears as a function of 
centrality, until « 40-50% most central Au-i-Au collisions show a FB correlation strength consistent 
with only short-range correlations. The evolution as a function of Arj in semi-peripheral data 
shows a similar behavior across all energies and systems, including pp. The energy dependence 
of pp shows that for energies greater than ^/s]^ = 200 GeV the FB correlation strength plateaus 
at a similar (small) value for Arj > 1.0, while y/sj^ = 62 A GeV data goes smoothly to b = 0. 
These long-range correlations can be ascribed to multiple elementary inelastic collisions, which 
are predicted in the Dual Parton Model and the Color Glass Condensate/Glasma phenomenology. 
The centrality dependence of the FB correlation strength may indicate the increasing dominance of 
a partonic core in the resulting particle production. 
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